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Abstract We employ ring-diagram analysis to study the sub-surface thermal 
structure of active regions. We present results using a large number of active 
regions over the course of Solar Cycle 23. We present both traditional inversions 
of ring-diagram frequency differences, with a total sample size of 264, and a 
statistical study using Principal Component Analysis. We confirm earlier results 
on smaller samples that sound speed and adiabatic index are changed below 
regions of strong magnetic field. We find that sound speed is decreased in the 
region between approximately r = 0.99 i?© and r = 0.995 i?Q (depths of 3 Mm 
to 7 Mm), and increased in the region between r ~ 0.97 i?© and r = O.985i?0 
(depths of 11 Mm to 21 Mm). The adiabatic index [Fi] is enhanced in the 
same deeper layers that sound-speed enhancement is seen. A weak decrease in 
adiabatic index is seen in the shallower layers in many active regions. We find 
that the magnitudes of these perturbations depend on the strength of the surface 
magnetic field, but we find a great deal of scatter in this relation, implying other 
factors may be relevant. 

Keywords: Helioseismology, Observations 
1. Introduction 

Measurements of the thermal structure beneath active regions using the tech- 
niques of local helioseismology are of substantial interest in solar physics. An 
increasingly wide variety of sophisticated models of the structure of sunspots 
and active regions are becoming available, and helioseismology can be used to 
test the validity of these models. In this work, we use ring diagrams to study a 
large number of active regio ns from So lar Cycle 23. 



The use of ring diagrams (jHilll . 119881) to study the near-surface acoustic prop- 



erties, dynamic s, and thermal st r uctur e of the Sun is by now well established 
(e.g. review by Gizon and BirchL 2005 ). In previous works, inversions of ring- 



diagram frequencies for structure have been performed on small numbers of 
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active regions ( Basu. Antia. and Bogart . 2004 : Bogart et al . 2008[ ) to determine 
the changes in sound speed and adiabatic index. In these works, sound speed and 
adiabatic index were found to be enhanced in the layers between approximately 
0.975 i?0 and 0.985 i?0, and depressed in the shallower layers between 0.99 i?© 
and 0.998 i?0. 

In this work, we study the sub-surface structure of a much larger sample of 
active regions than have prev iously been examined. We proce ed i n two ways: i) 



we sim ply extend the work of iBasu. Antia. and BogartI (j20Q4l) and lBogart et al. 



(|2008f ). performing structure inversions on ring-diagram data in the same manner 
as these earlier studies, and ii) in order to better take advantage of the large data 
set that we have, we perform a Principal Component Analysis (PCA) and study 
the resulting principal components. In Section [2] we introduce the ring-diagram 
data sets that we use in this work and discuss the differences between the sample 
that we use for individual inversions and the sample that we use in our PCA. 
In Section |3l we describe the use of Principal Component Analysis to study the 
ring-diagram mode parameters. In Sectional we present the inversions of both a 
large number of individual rings and of the principal components from Section 
[21 We summarize our findings in Section [S] 



2. Data 

2.1. Ring Diagrams 

Ring diagrams are three-dimensional power spectra of localized regions on the 
Sun. In this work, the data are resolved line-of-sight Doppler measurements 
( "Dopplcrgrams" ) taken with the Michelson Doppler Imager (MDI) instrument 
on the Solar and Heliospheric Observatory (SOHO). In order to achieve suffi- 
ciently high spatial frequency and sufficient signal-to-noise, we use only full-disk 
Dopplcrgrams, and require a duty cycle greater than 80%. Because full-disk data 
arc an MDI high-rate data product, sufficient data coverage is only available 
during yearly Dynamics campaigns, which arc typically two or three months 
long. The data have a one-minute cadence and are tracked for 8192 minutes. The 
midpoint of the tracking interval coincides with the region of interest crossing 
the central meridian of the solar disk. The data are projected onto a rectangu- 
lar grid using Postel's projection and corrected for t he distortion i n the MDI 
optics. This process is de s cribed in more detail by IPatron et al. ( 1997 ) and 
Basu. Antia. and Tripathv ( 1999() . 



The spectra are fit in the way described bv lBasu. Antia. and BogartI (|2004[ ) . A 
model of the power spectrum with 13 free parameters is fit at constant frequency. 
This model accounts for advection in the zonal and meridional directions as well 
as asymmetry in the spatial frequency direction and azimuthally around the 
ring. The parameter of interest in this work is the spatial frequency [k] , which is 
returned as a function of the temporal frequency [v] and the radial order [n] . In 
the plane-wave approximation the spatial frequency is related to the spherical 
harmonic degree [i] by ^ = fci?© . We interpolate the ring fits to integer values of 
i and interpret them in the same way as global mode parameters. We will refer 
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to these interpolated ring-diagram fits as "mode parameter sets" ; all results in 
this work are based on these data sets. 

The level of activity in a ring is ch aracterized by a Magnetic Activity Index 
(MAI: iBasu. Antia. and Boeart . 2004 ). This is a measure of the strong unsigned 



magnetic flux within the ring-diagram aperture, averaged over the tracking 
period. 

2.2. Ring Diagram Selection 

In this work, we use two different data sets to study the structure beneath active 
regions. The first set consists of a set of frequency differences between active re- 
gions from the NOAA catalog and nearby regions of quiet Sun, and is essentially 



equivalent to the data sets used in earlier works (e.g. iBasu. Antia. a.nd Tripathv , 
iRaiaguru. Basu. and Antial . l200l [ iBasu. Antia. and Bogartl . |2004 [). albeit 



with a much larger sample. The second set consists of mode frequencies for a set 
of largely non-overlapping rings, and is used for the statistical study described 
in Section [21 

The first sample is the same sample used by Baldner. Boeart. and Basu ( 2011ah .| 



Regions are selected from the NOAA active-region catalog between June 1996 
and April 2008. We include only NOAA active regions that arc identified at disk 
center — in other words, we require that the active region be present through 
the center of the t racking interval. We requ i re, to o, sufficient data coverage as 



noted above. As in lBasu. Antia. and Bogartl (|2004[ ). we perform inversions rela- 



tive to quiet-Sun ring-diagram mode parameters. The comparison ring diagrams 
are tracked at the same latitude as the active region that they are compared 
with. They are tracked across the disk in exactly the same way as the active 
regions, and are located within 60° longitude of their active region. This is done 
to minimize the systematic errors in mode-parameter estimation arising from 
projection effects and secular changes in the instrument, as well as to minimize 
the effects of difficulties in modeling the near-surface layers of the solar interior. A 



prelimi nary version of this analysis was presented bv lBaldner. Bogart. and Basu 



( 2011bl ). Although the parent active- region sample is the same in both works 



the actual ring-diagram data used in the inversions have changed somewhat — 
in particular, the mode sets used and the comparison regions chosen. The most 
significant changes to the data set are improved treatment of the errors in the 
mode fitting and inversion, and more careful comparison region selection. 

The second data set in this work is used in the PCA described in the next 
section. We use mode frequencies rather than frequency differences for the PCA 
itself. The PCA requires a filled matrix of mode parameters [i'n,e] for a set of 
modes common to all of the rings in the sample. A common set of [n, i] modes 
are chosen, and the ring-diagram fits are interpolated to these targets. 

For the PCA procedure, we wish to minimize overlapping data, which could 
lead to some spurious statistical results. We choose only ring diagrams that were 
constructed from tracked data patches that do not overlap. The ring diagrams in 
this sample are the same size and tracking length as in the data set used for our 
individual inversions but are distributed somewhat differently so as to minimize 
spatial overlap rather than guaranteeing coverage of all NOAA active regions. 
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3. Principal Component Analysis 

The use of Principal C omponent Analysis (PCA) in the analysis of helioseismic 
data was described bv iBaldner and Basul ( 2008h in the context of global- mode 



analysis. In brief, PCA finds an efficient representation of a set of data vectors 
as a linear combination of orthogonal principal components. It is efficient in the 
sense that, when properly ordered, each principal component is responsible for 
less variance in the data than the one that precedes it. In many applications, the 
data can be adequately represented using a relatively small number of principal 
components. An added benefit is that errors in measurements can be greatly 
reduced. In this work, we are primarily interested in using PCA to separate 
spurious signals in the mode parameters from those actually associated with 
solar activity, and to parametrize the changes in solar structure around active 
regions in a simple manner. 

We use the following notation: if the set of all m frequency measurements 
lyj for the ith ring diagram form the vector — {i/q ■ ■ ■ i^m}, each observation 
Dj can be reproduced completely with the m principal components £,k (each 
a vector of length m) and the m scaling coefficients Ck^i associated with that 
measurement: 

rn 

D^ = 5IcM?fc- (1) 

fc=0 

The principal components [^^J are normalized to unity. Errors in each component 
of the principal components and in the scaling coefficients are computed using 
a Monte Carlo simulation. 

The results of the PCA are shown in Figure [TJ The first component may 
be interpreted as the basic £-1^ diagram for a ring-diagram fit. The scaling 
coefficients [ci] give the appropriate scaling for each observation, and we can see 
a general trend with increasing magnetic activity, although the scatter is large. 
The second component also has a significant dependence on magnetic activity. 
The following two principal components which are shown in Figure [T] and are 
dominated by tails at the end of each n ridge. These tails arise in the fits as 
the signal-to-noise decreases at the high-^ end of a ridge or the ridges begin 
to overlap significantly at low £. In cither case, the fits become less reliable. In 
inversions on real data, we truncate the fits as these tails become significant. 

Higher-order components are not shown — they are in general dominated 
either by tails or by small numbers of outliers. Most of these components are 
significant only in one or at most a small number of rings. Some components show 
dependence on latitude (for example, ^5 and ^12) or on epoch (^7). By choosing 
to neglect these components when reconstructing the data sets, we remove many 
small systematic effects and errors from the frequency measurements. 

Given the apparent dependence of the scaling coefficients ci and C2 on mag- 
netic activity, it is desirable to reduce these coefficients to simple functions of 
MAI. In Figure [TJ we show a smoothed fit using B-splines to the ci coefficients 
as a function of MAI. These may be used to approximate the dependence of 
these coefficients on magnetic activity. Since ci and C2 both depend on MAI, we 
examine in Figure [2] their dependence on each other. It is clear that they are 
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Figure 1. The first four principal components and the associated scaling coefficients. The first 
four principal components (^i through §4) are shown in the left-hand panels as a function of 
degree [£] . Each order [n] is shown in a different color. The right-hand panels show the scaling 
coefficients for each ring as a function of the ring's MAI. The red line in the top-right panel is 
a smoothed fit to the scaling coefficients. 



tightly correlated. Thus, for a given ci, which wc can relate approximately to 
MAI, we can determine C2 fairly accurately. The scaling coefficients C3 for the 
third principal component do not correlate well with ci or C2 , or with any other 
quantity that we examined. 

It is important to know how well we can reproduce the actual ring-diagram 
mode parameters with a small number of vectors. In Figure |3] wc plot the 
residuals between the reconstruction using three principal components of three 
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Figure 2. The second principal component scaling coefficients as a function of the first. The 
points are color-coded by MAI. The magenta line is a smoothed B-spline fit to the points. 



randomly selected rings, and the actual ring data. We find that, especially for low 
MAIs, the agreement is fairly good. At the edges of ridges there are discrepancies, 
but as noted above, the spectrum fits tend to become unreliable at either end of 
the power ridges. 



4. Inversions for Structure 

To determine the structure of the solar layers beneath active regions (or, in 
general, regions of high surface magnetic field), we invert frequency differences 
for sound speed or for the first adiabatic index [Fi]. The treatment of this 
problem as a linear inversion is by now well established both in global and in 
local sei smology. In this work, we use Subtra ctive Optimally Localized Averages 
(SOLA: IPiipers and Thompson! . Il992[ Il994[) to perform the inversions. 



In our SOLA inversion, there are four free parameters: the error suppression 
parameter (/i), the cross-term trade-off parameter [/3], the width of the target 
kernel [A^i], and the number of B-splines used to remove the surface term 
[A]. The selection of these parameters has been describe d in earlier works, e.g. 
Rabello-Soares. Basu. and Christensen-Dalseaardl ( 19991 ). 



4.1. Inversions of Individual Rings 

The principal difficulty in performing an inversion for structure from helioseismic 
data is in c hoosing the appropriate inversion p arameters, which has limited 



aata is m c noosmg tne appropriate inversion p arameters, wmcn nas iimitea 
earlier works. Baldner. Boeart. and Basul ( 2011lJ ) found that a few combinations 



of inversion parameters worked well on a small subset of the active regions in our 
sample that was studied in detail. In this work, we again performed a thorough 
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Figure 3. The difference between frequencies obtained from three ring diagrams and their 
representations using the first two principal components. The residuals are divided by their 
errors. The MAIs of the three rings are shown in the lower right. As in Figure [T] different 
orders n are plotted in different colors. 

exploration of the parameter space on a small subset of the rings, including 
comparisons to rings published in earlier works. We find that the inversion 
results are not strongly dependent on the choice of surface term, so long as 
a surface term is used. When the surface term is not included in the inversions, 
very large perturbations are returned by the perturbations which are almost 
certainly spurious. We restrict ourselves to A = 4. In exploring the effects of the 
target kernel width, we find smaller values of A(ro) tend to cause oscillatory 
solutions to many inversions. Since we have not found any particularly sharp 
features in our inversions, we use a fairly large value of A(ro) = 0.055, which 
suppresses some oscillatory behavior in certain inversions, and does not overly 
smooth actual structure in better inversions. 

Finally, for the choices of the error suppression term and the trade-off param- 
eter, we have found that the values of these parameters for acceptable inversions 
in the rings we have examined carefully fall in a fairly narrow range. For the 
full sample, then, we run a batch of inversions for each region over the sample. 
Examining the inversions for each region can then be done fairly quickly and 
the best set of inversion parameters selected. We reject regions with unstable 
inversions — that is, regions whose inversion results are very strongly dependent 
on the choice of inversion parameters. What remains comprises our sample of 
sound speed and adiabatic index inversions. 



SOLA: ms.tex; 3 October 2012; 0:10; p. 7 



C.S. Baldner et al. 



0.02 - 
0.01 - 
0.00 



-0.01 - 



-0.02 - 

-0.03 
0.02 

0.01 

CM 

^ 0.00 

CN 

"^D -0.01 



-0.02 



-0.03 




AR9899 

/\MAI = 35.4^ 





AR9901 

/\MAI = }73.9G 



0.97 



0.98 



0.99 



0.97 



0.98 0.99 

r/Re 



Figure 4. Examples of inversions for adiabatic squared sound speed [c^] for four active regions. 
Tlie sense of the inversions is active minus quiet. The plotted horizontal error bars are the 
distance between the first and third quartile points of the inversion averaging kernels, and 
represent the resolution of the inversion. The vertical error bars are the formal errors in the 
inversions. 



Inversions for the difference in squared sound speed [c^] and adiabatic index 
[Fi] were performed for all regions in the sample with A MAI > 40 G. Figure 
U] shows example sound-speed inversions for four different rings with a range 
of active-region strengths as a function of depth. Figure [5] shows inversions for 
adiabatic index for the same regions. 

Figure IHl shows averages of the inverted sound speed for all regions in our 
sample for different depth ranges. For sound speed averaged between 0.975i?o 
and O.985i?0, sound speeds are generally enhanced in the presence of magnetic 
fields, while in the region from 0.99i?o to 0.998i?o, sound speeds decrease. In 
both regions, the magnitude of the change tends to increase with magnetic-field 
strength, although the relationship seems to be more of an envelope than a linear 
relation. Further, there seems to be some saturation of the effect at very high 
magnetic-field strengths. 

In the shallowest layers (r > 0.998 Rq), the sound-speed inversion results show 
a sharp positive change. This feature is present in every inversion in our sample. 
Caution should be used in interpreting this, however, as the averaging kernels 
become strongly asymmetric near the surface and the cross-term contributions 
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Figure 5. Examples of inversions for adiabatic index [Fi] for four active regions. The regions 
shown are the same as those in Figure 2] 
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Figure 6. Averages of inversions for over two depth ranges are shown, plotted as a function 
of A MAI. Blue points are averages of inverted sound speed between 0.975 Hq and 0.985 ii0; 
red crosses are averages of inverted sound-speed difference between 0.99 iJ© and 0.998 /?0. 
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Figure 7. Averages of inversions for differences in Fi over tv^o deptfi ranges are showrn, plotted 
as a function of A MAI. Blue points are averages of inverted sound speed between 0.975 -Rq 
and 0.985 Rq; red crosses are averages of inverted sound speed between 0.99 Rq and 0.998 Rq. 



begin to become significant. Examples of this sharp feature can be seen in the 
inversions in Figure SI 

We have also inverted for the first adiabatic index. In Figure [71 we show the 
same averages as Figure [SJ but for the Fi inversions. In general, Fi is an easier 
quantity to invert for (in the sense that the inversions tend to be less sensitive 
to the inversion parameters) , and so we include a larger number of inversions for 
Fi than we could for c^. 

We find that, as for the results, there is a depression in Fi in the shallower 
layers that we invert (above r = 0.99 i?©), and that, in many cases, there is a 
corresponding enhancement below approximately r = 0.98 i?©, as was found in 
earlier works. We find that the deeper enhancement is, for most rings, much less 
pronounced than for the c^. For some regions, in fact, we do not see any positive 
perturbation at all, and in general we find only a weak correlation with magnetic 
activity. 

The depth ranges of the negative and positive perturbations are found to 
be relatively constant. In Figure [SI we plot the depths at which our inversion 
results change from negative to positive. We find that, at lower activity levels, 
there is greater variance in the depths of the perturbations, but this can be 
explained by lower signal-to-noise in the inversions. Beyond that, we do not find 
any significant change in the depth ranges of the perturbations. 
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Figure 8. Position of transition, in fractional radius, between the negative (shallower) in- 
version results and the deeper (positive) perturbations in the inversion results, plotted as 
a function of A MAI. At left, the sound speed results are shown, at right are the results 
for adiabatic index. These transition points are found by treating the inversion points as a 
continuous curve and finding the zero-crossing point. 



4.2. Inversions of PCA Data 

In this section we use the PCA reconstruction of ring-diagram mode param- 
eters described in Section [3] to perform inversions for structure. As discussed 
earher, the advantages of using PCA on this data set are reduced errors in the 
reconstructed data set and removal of certain systematic effects in parameter 
estimation due to secular changes in the MDI instrument over time and projec- 
tion effects. In addition, because we may reduce the data to a number of linearly 
independent vectors, we substantially reduce the number of inversions that need 
to be done. 

Most of the variation in the ring-diagram mode parameter sets are spanned by 
the first two principal components (see Figure[3]). Furthermore, as the coefficients 
Ci are relatively well correlated with MAI, and C2 is tightly correlated with ci, we 
can reconstruct ring-diagram frequencies as a function of MAI. We then compute 
frequencies [t'n,^], the difference in frequency [i^t'n^f] between the target MAI and 
zero MAI, and the error in Svn,e- With these quantities, we may perform a SOLA 
inversion in the same manner as in the previous section. 

The inversion results for six target MAIs between 50 G and 300 G are shown 
in Figure m As in the inversions of individual rings, the inversions of the PCA 
reconstructions show a distinctive two-layer structure with a shallow negative 
perturbation and a deeper positive perturbation. In addition, the magnitudes of 
these perturbations scale with MAI, and the saturation effect hinted at in the 
results in Figure [5] is clearly present here. 

While most of the variation is spanned by and ^2 , there remains significant 
signal in subsequent components. We are interested to know, then, what effects 
these higher components have on inversions for structure, and so we invert 
components separately. It must be noted, however, that inversion results from 
different components, each with different errors, cannot be simply combined 
together linearly. With that in mind, we proceed to invert various principal 
components separately. 
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Figure 9. Inversion of the first two principal components at diff'erent MAIs for c^. The scaUng 
coefficients are chosen as discussed in the text. The left-hand panel shows individual inversions 
as a function of radius. As in Figure 2] the horizontal error bars are the difference between 
the first and third quartile points of the averaging kernels, and represent the resolution of the 
inversion. The vertical error bars are the formal errors in the inverted quantities. The latter 
are generally smaller than the points. The right-hand panel shows the average of the inversion 
results over the same depth ranges as Figure [6] as a function of MAI. 



In Figure [TUl we show inversions for sound speed for principal components ^2 
through ^5. Because it is not clear how to parametrize the scaling coefficients 
for any of the components except ^2, we simply scale the component by the 
largest coefficient for that component, and perform the inversion. We invert 
the components themselves, not explicit differences — thus the inverted sound 
speeds are the difference between a ring with a large coefficient Ci and a ring 
without any contribution from that component at all. We see that ^3 and ^4 do 
return changes in structure, but ^5 is consistent with no change. 

We also invert for adiabatic index. Figure [TT] shows inversions for the first two 
principal components at various MAIs between 50 G and 300 G. We find general 
consistency with the individual ring-diagram inversions; a shallow negative and 
deeper positive perturbation. The correlation with MAI in the shallow layer is 
very weak, and in some rings there is no negative perturbation at all. We find 
that the boundary between the negative and positive perturbations is somewhat 
deeper than in the individual inversion for Fi, occurring between 0.975 i?Q and 
0.98 i?©. This discrepancy does not appear to arise from the choice of modes 
used in the inversion, nor does it depend on the choice of inversion parameters. 
One possibility is that the contributions from other principal components may 
shift the boundary upward, which would in turn imply a dependence of Fi 
on something other than MAI. This is in fact what we find. On examination, 
the ^4 principal component returns a positive perturbation in the region we 
expect, and little signal elsewhere. The magnitude of the ^4 component (the 
C4 coefficients from Figure [T|) show a large scatter at low MAI values, but a 
pronounced bump above approximately 150G. This is broadly consistent with 
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Figure 10. Inversion of the §2 through §5 principal components individually. Each component 
5i is scaled by the largest associated coefficient Ci in our data set. The base frequencies are the 
same as the base frequencies used in Figure (9] 

the individual inversion results shown in Figure [71 where many of the regions 
with A MAI < 150G do not show positive perturbations. 

We have explored the first ten principal components in some detail and not 
found evidence of this, however. 



5. Conclusions 



In this work, we have presented inversions for structure of a large number of 
active regions, and used Principal Component Analysis to parametrize more 
accuratel y these inversion results. We confirm earlier results u sing ring-diagram 
analysis ( Basu. Antia. and Bogarti 20041 Bqgart et all . 12008 ). as well as high- 
degree global-mode analysis (jRabello-Soaresl 120121 ) that finds a two-layer ther- 
mal structure beneath surface magnetic activity. We find a shallow negative 
sound speed perturbation and a somewhat deeper positive perturbation, with a 
similar structure in the adiabatic index. The magnitudes of these perturbations 
ge nerally increa s e with the magnitude of the surface magnetic field. 



Bogart et al\ (j2008[ ) found a linear correlation between magnitude of the 



sound-speed change and strength of the active region. We find a similar relation, 
but with substantial scatter in our inversions of individual rings. It is perhaps 
more correct to say that we find an "envelope" related to MAI than a real 
correlation. Further, we find that the correlation appears to saturate at high field 
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Figure 11. Inversion of the first two principal components at different MAIs for adiabatic 
index [Fi]. The data are the same as those inverted in Figure |9] The left hand panel shows 
the individual inversions, the right hand panel shows averages over the same depth ranges as 
Figure [7] as a function of MAI. 



strengths. We find that the positive perturbation in Fi is much weaker than in c^, 
and that the correlation with magnetic activity is much less significant. We find 
consistent values for the s t rengt h of the pe rturbations in the sh allower layers with 
iBasu. Antia. and Bogart ( 2004 ) and with Bogart et al. ( 2008 ). but we find con- 
sistently smaller values for the deeper, positive perturbations than those found 
in the previous works. The choices of inversion para meters and mode sets can 
have s ome effect on the magnitudes of the inversions. In lBasu. Antia. and Bogart 
(j2004l ). RLS inversions tended to return slightly larger perturbations than the 
SOLA inversions, which we employ in this work. The most significant effect on 
our results appears to be our decision to neglect the /-mode in our inversions, 
which, as noted above, made our inversion results more stable. It also appears 
to have decreased the m agnitude of the sound-speed change in the deeper layers. 
Bogart et al. (j2008[ ) also reported that the depths of the positive perturba- 



tions in Fi were deeper than those in the inversions. Wc find this unambigu- 
ously in the PC A inversion results (Figure [TT|) . and it is also consistent with the 
inversions of individual rings. 

Interpreting these results is not straightforward. The inversions that we per- 
form in this work are, strictly speaking, only valid in a spherically symmetric, 
non-magnetized star. Clearly, neither of these conditions is satisfied in active 
regions. In the presence of magnetic fields, wave propagation is affected both 
by the direct effect of the Lorentz force, and by the effects of magnetic fields 
on the thermal structure. At some point — at some minimum field strength or 
field configuration — the entire linearized inversion becomes meaningless, as the 
contributions from magnetic fields to the wave propagation (which is, among 
other problems, no longer insensitive to direction) become significant. 
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Structure of Active Regions 



In some studies, the choice is made to interpret the infer red change in sound 



spe ed as a change in th e local wave speed, as was done in Gizon et al. ( 20091 ) 



and lMoradi et all (|2010l ). rather than interpreting it as the effects of a thermal 



perturbation as we do here. In these works, significant disagreement was found 
when ring-diagram inversions were compared to time-distance inversions for the 
same solar data. The resolution of this discrepancy re mains unclear. If w e do in 
fact measure a wave speed perturbation in this work, Lin. Basu. and Li ( 20091 ) 
have claimed that the magnetic effect can b e disentangle d from the thermal 
effect. In the case of time-distance analysis, Braun et al. ( 2012[ ) have shown 



that travel-time shifts determined from a model sunspot match neither what 
the thermal perturbation should give nor what they would expect from the 
magnetoacoustic fast mode speed. The implication of this work is that, in the 
case of time-distance analysis, at least, it is not correct to interpret inversion 
results as either changes in the thermal sound speed or as local wave speed 
perturbations. Equivalent work has not yet been done for ring-diagram analysis, 
but must be done to determine the extent to which the assumptions we have 
made are valid. 

In studying our sample of rings using PGA, we have decomposed the ring 
diagram frequency measurements into a set of linearly independent components. 
The dependence of the first two of these on magnetic activity allows us to 
parametrize most of the frequency variance across our sample of rings as a 
function of MAI, and to determine what changes in sound speed and adiabatic 
index could give rise to these changes. We find that both sound speed and 
adiabatic index have a dependence on MAI that is consistent with what has 
been found in individual rings. 

There is further variance in the ring-diagram frequency measurements, how- 
ever, and MAI alone does not adequately parametrize the ring-diagram frequen- 
cies, as can be seen in Figure [TJ Many of these changes are due to errors in the 
ring fitting and systematic effects due to projection effects and secular changes 
in the MDI instrument, but some of the more significant principal components 
may represent changes that are solar in origin and that may be associated with 
thermal changes below the solar surface, as shown in Figure fTOl 

We have demonstrated that PGA can be a useful tool in ring-diagram anal- 
ysis and structure inversion. The large volume of data being returned from the 
Helioseismic and Magnetic Imager (HMI) on the Solar Dynamics Observatory 
(SDO) spacecraft represents a significant data analysis challenge — it is possible 
that this technique might prove a feasible alternative to attempting inversions 
on every ring diagram produced by HMI. 
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